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Abstract. Corrections induced by primordial non-Gaussianity to the linear halo bias can 
be computed from a peak-background split or the widespread local bias model. However, nu- 
merical simulations clearly support the prediction of the former, in which the non-Gaussian 
amplitude is proportional to the linear halo bias. To understand better the reasons behind 
the failure of standard Lagrangian local bias, in which the halo overdensity is a function of 
the local mass overdensity only, we explore the effect of a primordial bispectrum on the 2- 
point correlation of discrete density peaks. We show that the effective local bias expansion to 
peak clustering vastly simplifies the calculation. We generalize this approach to excursion set 
peaks and demonstrate that the resulting non-Gaussian amplitude, which is a weighted sum 
of quadratic bias factors, precisely agrees with the peak-background split expectation, which 
is a logarithmic derivative of the halo mass function with respect to the normalisation am- 
plitude. We point out that statistics of thresholded regions can be computed using the same 
formalism. Our results suggest that halo clustering statistics can be modelled consistently 
(in the sense that the Gaussian and non-Gaussian bias factors agree with peak-background 
split expectations) from a Lagrangian bias relation only if the latter is specified as a set 
of constraints imposed on the linear density field. This is clearly not the case of standard 
Lagrangian local bias. Therefore, one is led to consider additional variables beyond the local 
mass overdensity. 
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1 Introduction 

Recently, galaxy clustering has become a serious contender to the cosmic microwave back- 
ground (CMB) in the hunt for primordial non-Gaussianity (see [1, 2] for recent reviews). 
Non-Gaussian initial conditions can indeed leave strong scale-dependent signatures in the 
galaxy power spectrum [3] and bispectrum [4]. However, one of the main caveats in the the- 
oretical modelling of galaxy clustering statistics and the interpretation of large scale surveys 
data is galaxy biasing: unlike CMB temperature measurements, galaxies furnish a distorted 
picture of the primordial curvature perturbations because they preferentially trace overdense 
regions of the Universe. This induces a bias between the galaxy and matter distributions, as 
was first pointed out in [5]. 

One of the simplest biasing schemes is the local bias model [6]. In its Lagrangian for- 
mulation, galaxies are a Poisson sampling of some continuous overdensity field dg{x) which 
is a function of the local, linear matter overdensity d{x) only. In general, 5g{x) is writ- 
ten as 6g{x) = bi6{x) + b26'^{x)/2 + • • • , where the bi are the Lagrangian bias parameters. 



- 1 - 



For Gaussian initial conditions, the predictions of this model are in reasonable (though not 
perfect) agreement with halo clustering statistics extracted from N-body simulations once 
gravitational evolution is taken into account. However, for a non-zero bispectrum of pri- 
mordial curvature perturbations, the local bias expansion gives rise to a correction to the 
linear halo bias whose amplitude is proportional to 62/NL [7~9], whereas N-body simulations 
clearly indicate that the ampltitude scales according to &i/nl [3, 10, 11]. By constrast, the 
peak-background split ansatz [5] predicts the correct amplitude oc bi for the non-Gaussian 
bias [12, 13], as well as an additional, initially overlooked correction [14, 15] that significantly 
improves the agreement with N-body simulations [16]. Does this imply that we should give 
up on local bias when it comes to primordial non-Gaussianity? Hopefully not, especially 
since a working local bias approach would be extremely useful for the computation of the 
galaxy bispectrum with generic initial conditions [8, 9]. 

Several alternatives to local Lagrangian bias have been explored in the context of pri- 
mordial non-Gaussianity. In the approach of [17], which is an extension of [18], the halo 
overdensity field is a function of both the local matter density and (non-local) Gaussian 
part of the primordial Bardeen potential ^(x). However, the non-Gaussian bias acquires 
additional corrections that strongly depend on the smoothing scale Ri (i.e. the cell size, 
which is distinct from the Lagrangian radius Rs of a halo). To cure this problem, [19] have 
recently proposed a renormalization procedure that absorbs i?;-dependent terms, so that 
the non-Gaussian bias amplitude is effectively independent of Ri. Furthermore, statistics of 
thresholded regions can provide insights into the relationship between local bias and peak- 
background split [20] . In the specific case of sharp clipping, both peak-background split and 
the corresponding local series expansion consistently predict the non-Gaussian halo bias [21]. 
However, one of the main drawbacks of this model is its poor fit to the clustering of dark 
matter halos with Gaussian initial conditions. Clearly, it is essential that the model be able 
to reproduce Gaussian clustering statistics as well. 

In this work, we compute the non-Gaussian correction to the linear bias within the 
peak formalism [22]. The central assumption of this model is the correspondence between 
virialized halos and maxima of the initial density field, so that the biased tracers constitute a 
genuine point set. Another essential difference with most of the approaches mentioned above 
(thresholded regions excepted) lies in the fact that clustering statistics of initial density 
peaks are fully specificed upon enforcing the peak constraint, without the need for local 
bias expansion or peak-background split. Yet another motivation for our work is the non- 
universality of the multiplicity function of initial density peaks. Non-universality of the 
mass function must be accounted for to obtain exquisite agreement with measurements from 
N-body simulations [15]. 

The outline of the paper is as follows. First, we demonstrate that the non-Gaussian 
bias obtained from a direct calculation of the non-Gaussian 2-point probability density agrees 
with that derived from the effective, local bias expansion advocated by [23]. This shows 
that the latter can be extended to non-Gaussian initial conditions. Second, we prove that 
the resulting non-Gaussian bias agrees with that inferred from peak-background split, in 
which primordial non-Gaussianity rescales the amplitude of density fluctuations erg in a scale- 
dependent manner. Third, we generalize these results to excursion set peaks [24] and, in 
particular, provide a local bias expansion along the lines of [23] that includes the first- 
crossing condition. We find again consistency with peak-background split expectations. This 
is especially interesting since excursion set peaks have recently been shown to give a very 
good match to the Gaussian mass function and linear bias of dark matter halos [25] . Finally, 
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we discuss our results in light of recent work before summarizing our findings. 



2 A brief overview of the peak formalism 

We are interested in the 2-point correlation function of density maxima of the initial density 
field S{x), which may not follow perfect Gaussian statistics. The assumption is that nonlinear 
structures form from the gravitational collapse of local maxima of the linear, smoothed 
overdensity field ^^(a;). To define a local maximum, one must require that the first derivative 
vanish and the second derivatives be negative. Following [22], we normalize 6s{x) and its 
derivatives as follows: 

u{x) = —6s{x) , r]i{x) = —di5s{x) , Cij{x) = —didj6s{x) . (2.1) 

O-Q O"! CJ2 

Here, the spectral moments aj are given by 

d^k 1 



where Ps{k) = W'^{kRs)P5{k) is the power spectrum of the initial density field smoothed 
on the Lagrangian scale Rs of a halo with a spherically symmetric kernel W{kRs). Density 
maxima are a point process and, therefore, their abundance is formally a sum of Dirac 
distributions. On using the Kac-Rice formula [26, 27] and enforcing the peak constraint, the 
number density of peaks of height Vc = Sc/ctq, where 5c is the critical threshold for spherical 
collapse [28], in a 3-dimensional Gaussian random field can be written as [22] 

33/2 

npk{x; Vc) = IdetCI Soin) Onih) 5d{v - i^c) ■ (2.3) 

where i?* = \/3o"i/cr2 is the characteristic radius of peaks and A3 is the smallest eigenvalue of 
the matrix —Cij, so that the term 9h{^3) ensures that we are only counting density maxima. 
Here and henceforth, we will omit the dependence of npk(x; Vc) on Rg for concisencess, but 
one should bear in mind that depends explicitly on both z^c and Rg. Note also that 
dealing with peaks of the density field does not imply V$ = at the peak positions. This 
assumption was sometimes invoked in the derivation of the non-Gaussian bias [3, 17] 

The n-point correlation function of peaks of given significance Uc is defined as the en- 
semble average 



-1 ^ _ {np]^{xi;uc)np]^{x2; Uc) ■ ■ ■ npk(a;„; Uc)) , , 

i+l;^^[XuX2,---Xr„Uc)- Mx; U^T ' ^ ' 



where the average is taken with respect to the derivatives of 6s{x), 

JJnpk(a^,;z^c)) = (HldetQ SdM OH[h(i)] Siit^i - i^c)) (2.5) 




3 ' \ f 



Here, Pn{xi;vi,r)iXi) is the joint probability distribution for the variables (z^, 7?,C) at n dif- 
ferent spatial locations. The n-point connected or irreducible correction function is defined 
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in such a way that {np]^(xi; v^) ■ ■ ■ n.pk(a^; Vc)) is a sum of terms, each of which pertains to a 
partition of the set of n locations xi, - ■ ■ (e.g. [29, 30]). In this paper, we will focus on the 
connected 2-point correlation of peaks of height Uc (at a separation r = \x2 — xi\) 

4?(r; .e) ^ ^-P^(-^^->P^(-^^^--)) - 1 . (2.6) 

This 2-point correlation function can be evaluated in two different ways, either from an 
explicit computation of {np\^{xi]Vc)np\^{x2;Vc)) for Gaussian initial conditions, see [22, 31- 
33]) or from the much simpler local bias approach proposed by [23]. Reassuringly, the two 
methods agree, yet the second is much faster and physically more intuitive. 



3 Peak 2-point correlation with non-Gaussian initial conditions 

In this Section, we derive the non-Gaussian corrections to the 2-point correlation function 
,^pj^ (r; fc) in two different ways: from a direct computation of the joint probability density 
P2(xi; i^i,r)i,Ci) and from the effective local bias approach [23]. We describe the two methods 
and demonstrate that they lead to the same non-Gaussian bias correction. 

3.1 Prom an Edgeworth expansion of the probabihty distributions 

In the presence of non-Gaussian initial conditions the generating function formalism can be 
applied to find an explicit expression for the joint probability Pn{xi; i/i,riiXi) before comput- 
ing correlators of npk(a;, z^c)- This functional approach has been used in, e.g. [34-37]. Here, 
we briefly review the basic arguments in the form appropriate for our purpose. 

3.1.1 Edgev^orth formahsm 

Let P{yi,y2r ' ' iUn) be a generic A^-dimensional probability distribution function. The 
connected correlation functions of the random variables yi can be written in terms of the 
generating function Z{Ji, J2, • • • Jn) = Z{J), which is defined as 



Z{J)= / dyidy2---dyNP{yi,y2,---yN)e 



jJlJ/l+iJ2J/2- 



-iJNVN 



{e'^-yi) 



(3.1) 



The connected correlators are then given by derivatives of the logarithm of Z{J), 

^ ^mi^ hmjv 



,(c) 



5mim2-m,jv 



/ mi n 



,m2 



df^'^djp ■■■dJ 



N 



logZ(J) 



Jl = J2 = 



=J«=0 



(3.2) 



Note that the yi may be defined at different spatial locations, but this is not essential to our 
discussion. The joint probability P{yi, ■ ■ ■ , yn) can then be constructed from the generating 
function Z{J). For random variables with zero expectations values, we have 



P{yi,y2, 



,yN) 



exp 



E 



.n=3 mi,---m]v 
1 



Ac) 

<,mi---mjv 

mil • • • vtinI 



dyi. 



d 



dyN 



V(27r)^detM 



exp 



(3.3) 
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where rrii = N and Mij = {yiyj) is the covariance matrix. We can replace the derivatives 
with respect to yi by introducing multi-variate Hermite polynomials defined as 

dyi) \dyN 



Hm{y;M~') = (-l)"exp Q/M^iy^ 



exp 



/M-iy , (3.4) 



where y = (yi, • • • , yjy) and m = {mi,m2, ■ ■ ■ niN)- On expanding the exponential factor in 
the correlations £,m\---mN with n > 3 (assuming they are small relative to the second order 
moments), we can eventually write down the joint probability density as 



P{yi,y2 ■■■ ,yN) 



1 



; exp 



V(27r)^detM 

1 + E E M y;M-') + 



n=3 mi,---mjv 



(3.5) 



The first line on the right-hand side is a multivariate Gaussian, whereas the second line 
is the non-Gaussian correction represented as an Edgeworth series in the n > 3 connected 
correlation functions. This series expansion can then be substituted in place of P„ in (2.5) 
to evaluate the non-Gaussian contributions to the A^-point function. 

3.1.2 Non- Gaussian bias 

We restrict ourselves to the corrections from the non-Gaussian sector of our interest, i.e. 
the contributions from 3-point correlation functions including the derivatives of the density 
field. A total of eight different combinations are possible, not including the exchange of the 
coordinates xi X2. Namely, these are 



{u^{xi)u{X2)) 
{u{xi)Cij{xi)u{X2)) 
{7]i{xi)r]j{xi)u{X2)) 
rn 



{v'^{xi)Cij{X2)) , 
{u{xi)Cij{xi)Cira{x2)) , 
{Viixi)Vjixi)Clm{X2)) , 
{Cij{xi)Ckl{xi)Cmn{x2)) 



(3.6) 



These 3-point correlators can generically be written as {Ti5s{xi)T25s{xi)T35s{x2)) , where Tj 
are differential operators. For instance, setting Tj = — yields the curvature u = —ASs/cr2 
of the smoothed density field. After some algebra, the Fourier transform of this 3-point 
correlator is 



d'r {Ti6,{xi)T25s{xi)T35,{x2)) e 



-ikr 



(2vr)= 



T^{q)T2{-k- q)T,{k)Ms{q)Ms{\k+ q\)Ms{k)i'i\q,-k- q,k) . (3.7) 



Here, the 3-point function ^^^^(^1, ^2, ^s) is that of ^{x), and M.s{k) = M.{k)W {kRg) where 
M[k) is the transfer function between ^{x) and the linear, smoothed density field at redshift 
z. Here and henceforth, we will consider the limit A; — )• 0, in which the right-hand side is 
approximately 



T^{k)Ms{k)P^{k) 



d^q \KA-i( \S!iSl~3l~R 
-^Ti{q)T2{-q)M,{q) 



(27r)2 



(3.8) 
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Clearly, this expression scales proportionally to only if T^{k) = 1 or, equivalently, if 
T'i&si'^) = i'{x2)- In the specific case of the constant-/NL</'^ model, the operators Ti{±q) 
give rise to the spectral index or o"!, depending on their detailed form. In the general 
case, the scale-dependence and amplitude of the non-Gaussian bias correction is governed by 
the squeezed limit of the bispectrum of curvature perturbations. Accounting for a possible 
scale-dependence of /nl, we parametrize the latter as [14, 19] 



-k) 2AMkp) 



2ai 



2q2 



k J 



(3.9) 



Here, kp is some reference wavenumber. This parametrization encompasses many bispectrum 
shapes considered in the literature. For instance, the quadratic coupling /nl^^ with k- 
dependent /nlI^) oc k'^f considered in [38] (see also [39]) corresponds to A = 1, ai = nj /2 
and 02 = 0, whereas the equaliteral bispectrum shape [40] yields A = 3/2, qi = and 
Oi2 = {^s — 4)/3. Consequently, (3.8) simplifies to 



T^{k)Ms{k)P^{k)2A 



-k- 



(2^ 



q^(-^+-^)T^{q)U-q)MUq)PM , (3.10) 



and is proportional to the square fj? of a spectral moment (2.2) with j determined by the 

values of ai, Q2 and the scaling of the Fourier transforms Ti{q) and ^2(9). 

(2) 

The calculation of the non-Gaussian corrections to Cp]l{^]^c) from the Edgeworth ex- 
pansion is straightforward. Including symmetric numerical factors and denoting a = ai + a2, 
the contributions of the 3-point correlators given in (3.6) to the 2-point power spectrum of 
discrete density peaks are, formatted as in (3.6), 



b2o{l)hoi2)al , 62o(l)feoi(l)fcV, 



2^2 

a 1 



26n (1)610 (2)a2 
-36io(2)ct; 



+1, 26n(l)6oi(2)A:V2^i, 
a+i, - 36oi(2)A;^cr^+i , 



(3.11) 



and 



^021 



l)-5-29„logGi"\7i,7ii/ 
602(1) -5-2a„ log G;,"\7i, 711^) 



0=1 



610(2)0-^+2 > 



Q = l 



6oi(2)A:V; 



2,2 

a+2 ) 



(3.12) 



where we have omitted the common multiplicative factor 4^/nl(A:p)A:-2°i A;-2"2A^,(A;)P^(A:) 
for brevity. Explicit expressions for the linear and quadratic bias factors bij are given by 
[23, 32, 33] 



6 



10 



501 



520 



6?, 



-0^(1 



7?) 



^11 



b^bu + 



7? 



7?) 



502 



62 



(3.13) 



Here, 71 = o"2/(c7oO"2) is a dimensionless number that takes values between zero and unity 
depending on the shape of the smoothed density power spectrum [22]. Moreover, 



where the overhne designates the average over u at the locations of density peaks. The 
function g1"\7i, 71 i^c) is defined in [33], 

(n - zf 



f{u,a) 




exp 



dw + 



2(1-7)^ 



u/2 



dv 



u/4: 



3v—u 



F{u, V, w) = {u — 2w) [{u + w)'^ — 9v'^] v ^ — w^) 



dwj F(n,^;,u;)e-^"(3^'+'"')/^ 

(3.15) 



Furthermore, although the bias parameters bij do not depend on position, we have added an 
argument of (1) or (2) to indicate whether they were evaluated from variables at position xi 
or 032. Adding up all the terms, we eventually obtain 



APpk(A;) = AA ^-^^ k-^^'Ms{k)P^{k) {bio + hoik") 

X <; alb2o + 2cj2+i6ii + ^2^2^02 - 3ct2^i - [5 + 2 a« log g[,"\7i, 7iz^) 



(3.16) 







a=l. 







We recognize the linear peak bias ci{k) = {biQ+boik'^)W (kRs) and, inside the square brackets, 
the quadratic bias factors xio and xoi defined by [23], 



Xio 



3 

'2^' 



Xoi 



2al 



l + -9„logG("\7i 



a=l 



(3.17) 



We can read off the (scale-dependent) non-Gaussian correction Aci{k) to the linear peak bias 
from APpk(A;) = 2ci{k)Aci{k)Ps{k). Notice that, in order to highlight the connection with 
the formalism of Matsubara [20], we have adopted his CTv-notation for the renormalized bias 
functions. However, we emphasize that our calculation does not involve any renormalization 
scheme whatsoever. 

3.2 From a local bias approach to peak clustering 
3.2.1 Effective local bias expansion 

As shown in [23], the 2-point correlation function ^pj^ (r; Uc) of peaks of a Gaussian random 
field can be thought of, up to second order at least, as arising from the effective local bias 
expansion 

(^pk(a;) = aQhiQv{x) + a2boiu{x) + ^aQb2oi''^ (x) + aoa2bui/{x)u{x) + ^cr| 602 1*^(^3) 



(3.18) 



provided that we ignore all the contributions involving powers of zero lag moments. Here, 
r]'^{x) and C'^{x) are the square length of the vector of first derivatives and the trace of the 
squared traceless part of the hessian matrix didjSs{x) respectively, 



1 



r^\x)^-[V6,{x)]\ 



"2 



didj - ^5ijA)6s{x) 



(3.19) 



This expansion is local except for the filtering of the mass density field. It is effective in the 
sense that 5pk(a^) is not a count-in-cell overabundance of peaks, but merely some idealized 
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continuous overdensity field that can be used to derive Cpk ('^i ^c) without liaving to compute 
the joint probability density P2 that describes covariances between variables defined at two 
different locations. 

The peak bias factors bij and Xij can be derived from a peak-background split argu- 
ment where the long- wavelength perturbation shifts the mean of the probability distributions 
characterizing the rotational invariants ^{x), u{x), •rf'{x) and C^(^)- The variables v and u 
are distributed normally, whereas 3ry^ and 5C^ are distributed as chi-squared (x^) variates 
with 3 and 5 degrees of freedom, respectively. The resulting non-central probability densities 
can be expanded in sets of orthogonal polynomials which, on enforcing the peak constraints, 
yield the peak bias parameters. The orthogonal polynomials associated with the correlated 
variables v and u are bivariate Hermite polynomials Hij, whereas those pertaining to if'ix) 
and C^(^) are generalized Laguerre polynomials L^^^ . We thus have 



The probability density Pi can be factorized into the product of a bivariate Gaussian A/'(i^, u) 
times x|(3r/^) and xK^C^); where x^{x) is a ^^-distribution with k degrees of freedom [see 
also (5.6)]. 

The bias factors b^ can be easily computed from the generating function given in [23]. 
In particular, the linear and quadratic bias factors associated to the scalars v and u are given 
by (3.1.2), whereas those pertaining to the quadratic invariants are given by (3.17). Note 
that the term XoiC^(^) in effective local series expansion has a simple physical meaning: 
it represents the leading-order contribution of asphericity in the peak density profiles to the 
clustering of density peaks. Therefore, it vanishes in the high peak limit ^ 1 since high 
peaks are close to spherical. 

3.2.2 Non- Gaussian bias 

From the effective local bias expansion (3.18), the leading order contribution to the non- 
Gaussian peak 2-point correlation is immediately computed as 



The last-line in the right-hand side signifies that the other terms can be obtained upon replac- 
ing all occurences of ((Tq, 610, ^2) by (c2, &0I) '"2)- In the low-Zc limit, the various contributions 




(3.20) 



A^pk (^;^c) = 2fT||6io62o(i^ii^2) + 4fTocr26iofoii(z^iniz^2) + 2o-oo-|6io&02(^i?i^2) 

+ 4<To<TiXloblo(??ll^2) + 4croCr2XOlfelo(ClJ^2) 




(3.21) 



8 



to the non-Gaussian peak power spectrum are given by 

2 



d r{u {x2)y{xi)) e 



ik-r 



A 



hhikp 



k 



(3.22) 



fc-5.0 4 / fj, 



d3r(??2(a;2)z.(a;i))e-*^ 

d'r{e{x2Mx,))e-^''- ^(^YA^^^k-''^^Ms{k)P,{k) 



O"0 V 0-2 



The 5 correlators that involve u{xi) instead of I'ixi) yield an extra multiplicative factor of 
k"^. On summing all contributions, the leading-order non-Gaussian correction to the peak 
power spectrum reads 



APpk(A;) = AA ^^^^ k~^^-Ms{k)P^{k) {ho + boik^) 



(3.23) 



This result agrees exactly with (3.16) obtained from the Edgeworth expansion of the joint 
probability distribution P2{xi,] I'ijrfiXi)- We shall now establish a connection with the peak- 
background split ansatz. 



4 A peak-background split interpretation of the results 
4.1 Peak-background split and non-Gaussian bias 

As shown in, e.g. [12-15, 41], non-Gaussian corrections to the bias parameters can be com- 
puted using a peak-background split [5, 22]. In its simplest expression, the non-Gaussian 
contributions can be read off from a Taylor expansion of the halo mass function. For exam- 
ple, assuming that the Gaussian mass function of the tracers takes the form u^ScCTq) and 
considering the const ant- /nl*/"^ model for simplicity, the change induced by a long-wavelength 
perturbation {5i,<l)i) is 

dfi dft 

n [6c - Si,aoil + 2/nl0/), ^3, • • • ] n{5c, (Tq) - —61 + 2/NLa-o^— ^^'i H (4.1) 

oOc oao 

Therefore, the overabundance of tracers is given by 

« -^^6i + 2/nl (^P-) <Pi + ---^ h,5i + 2/NL6r<^/ + • • ■ (4.2) 

The non-Gaussian contribution to the first-order bias thus is proportional to the first deriva- 
tive of average number density n of halos with respect to ctq. However, while 6^*^ = 61 for 
a universal mass function, this will not hold in general. Is this simple intuitive picture also 
valid for the discrete peaks considered here? 
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4.2 Application to discrete density peaks 

To answer this question, we shall first apply the peak-background split to the peak number 
density npk, which is the ensemble average of (2.3). This was computed in [22] as 



1 



n 



pk 



(2^)2i?3 



(4.3) 



Interestingly, ?ipk is not a universal function as it depends distinctly on Sfi- and Rg- fact, 
we can also assume that it depends on ctq, cji and (T2 through the parameters 71 and R^, the 
significance fc and the normalized field u (which is integrated over in gE,°^). Hence, we can 
write rlpk = ^pk((5c, {cTj}), with i = 0,1,2. Therefore, the above discussion suggests that we 
take derivatives of npk with respect to the spectral moments ai. We begin with the derivative 
of npk with respect to ctq and easily obtain 



'^"-pk _ dup], duc ^ (9npk 571 
dao 



+ 



(4.4) 



Similarly, the derivative of ripk relative to the spectral moment ai is 
9npk _ 9npk dR^, 5npk 571 



dai dRi, dai 
Finally, the derivative of flpk with respect to cj2 yields 



= 2o-inpk(^Xio + oiij • 



(4.5) 



dn 



pk 



dao 



1 



du 







/(«) 



1 



+ 



dfju) 
dao 



-(«-7i^c)V[2(l-7?)] 



27r (1 



7?) 



J__ ^ dnp]^dR^ ^ dup^d-fi 
a2^^^ dRi, da2 d^i da2 



(4.6) 



Here, = (27r)^/^i?2 and the term — npk/c2 arises upon taking the derivative of the measure 
du with respect to a2- The derivative of f{u) = f{u,a = 1) produces a factor of — 8npk/cr2, 
to which the measure dvdw contributes — 2npk/c'"2 and the shape factor F{u,v,w) (which is 
a homogeneous function of degree —6 in CJ2) the remaining — 6npk/o"2. There is an additional 
contribution that originates from the dependence of 3f ^ + w'^ on a2 in the argument of the 
exponential. It can be written down as —{2/a2)daf{u,a)\a=i- Adding up all the terms, we 
find 



dn 



pk 



dao 



1 j — ^iVcU \ - 9 _ 

2 'bk "-pk 

0-2 V 1 - 7l / 0-2 



— a^log g[,"\7i, 71 i/c) 

0"2 



n 



a=l 



pk 



+ 



9npk dRi, dn^y, dj 



dRi, da2 

The logarithmic derivatives (npk/o-i)9npk/(?cjj thus are 

d log npk 



+ ^ — = cr2f^pk 2x01 
d-)i da2 V 



'02 



(4.7) 



^log^pk 2, , . 
<9 log ao 



d log cJi 



2(jf xio + 2a^ Oil , 



d log npk 

9 log (72 



20-2X01 +0-261 



'02 • 



(4.8) 
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On multiplying the logarithmic derivatives by (T^_)_j/cJ? and summing the resulting contribu- 
tions, we arrive at 



2 _ . -2 

y-v Olognpk / (Ta+i 



j=0 



d log Oi 



aih2Q + 2<+i6ii + ^^+2^*02 + 2(J^+iXio + 2<+2Xoi + ^ • (4.9) 



This is precisely the amplitude of the non-Gaussian bias correction found from the compu- 
tation of the peak 2-point correlation, except for an additional factor of ct^/iTq. 

4.3 From peaks to dark matter halos 

To explain the origin of this factor, we note that the differential number density of dark 
matter halos of mass M per unit comoving volume is generically expressed as 



nh(M) 



P ,1 ^ dloguc 



(4.10) 



where the multiplicity function /(I'c) encodes information about halo biasing. In the present 
calculation, the quantity V^npk(<5c, Wi}), where V = M/p is the Lagrangian volume associated 
with the filter (M oc Rg), plays the role of a multiplicity function. Therefore, we can write 
the halo mass function associated with the peak number density (2.3) as 



-h(M) = -npk(5e,{a.})^^ 



(4.11) 



The logarithmic derivative of nh(M) with respect to ai will be identical to that of npk except 
that, for (To, there will be an additional factor of —1 owing to the multiplicative factor of i^c 
in the expression of the halo mass function. This factor precisely cancels the extra term of 
(7^/iTo in (4.9). Therefore, we find 



2 _ . -2 

91ognh f aa+i 



i=0 



d log ai \ a, 



'2: '2i '2: 2 2 



(4.12) 



which is exactly the amplitude obtained from the calculation of the peak 2-point correlation. 
The non-Gaussian contribution to the linear peak bias thus is 



Acm = 2A^ 



2 — / \ 2 

d log rih ( (Ja+i 



E 

.j=0 



OlogfJi V a. 



(4.13) 



The physical interpretation is straightforward: in the presence of a primordial 3-point func- 
tion, a long-wavelength background perturbation of wavenumber k rescales the amplitude of 
the power spectrum Psiq) of the smoothed density field in a scale-dependent manner [12-14], 



Ps{q) 



1 + 2^)P.(,) 

0"8 



l + 2e 



kr 



2qi 



2a2 



Psiq) 



(4.14) 



where e will be determined shortly. Note that this approximation is valid in the limit k <^ q 
solely. In the particular case of the constant-/NL0^ model, all the spectral moments are 
rescaled analogously, iTj — )■ (1 + e)(Tj, so that the parameters 71 and R-^ remain unchanged. 
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In general however, the transformation ai ai + 6ai will be scale-dependent. We can write 
the derivative of with respect to the normalisation amplitude as 

^^Sas = y ( ^1 = efc-^- y f . (4.15) 

On setting e = 2f^i^{kp)Jv[~^{k), we recover the full non-Gaussian fc-dependent correction to 
the linear halo bias. 

5 Extension to excursion set peaks 

The differential peak number density ripk cannot really be interpreted as a multiplicity func- 
tion because it is defined for a fixed smoothing scale Rg, whereas one should allow Rg to vary 
while (5 = (5c is kept fixed. This is the reason why we have not yet recovered the strong mass- 
dependent correction found by [14] . Furthermore, the linear bias factor feio is not equal to the 
Gaussian peak-background split bias —d\n.n\i/ d6c inferred from a uniform shift (5c — ?• 5c + ei 
applied to the halo mass function (4.11). 

Therefore, this suggests that we consider the improved model of [24, 42], in which peaks 
on a given smoothing scale contribute to the multiplicity function only if the conditions 
5{Rs) > 5c and 5{Rs + Ai?s) < 5c are satisfied. In other words, the density must upcross the 
threshold for collapse on the smoothing scale Rs- Because the trajectory described by the 
Gaussian or tophat filtered 5{Rs) as a function of Rs is strongly correlated when (Tq{Rs) ^ 1, 
this almost certainly implies that 5[Rs) upcrosses the threshold for the first time at R = Rs 
(see, e.g. [43] for a detailed discussion). 

5.1 Effective local bias expansion 

To include the first crossing condition into the calculation of the non-Gaussian bias, we must 
introduce a new variable: /x = —d5s/dRs = —5's (this notation is borrowed from [42] who 
considered the particular case of a Gaussian filter). The number density of density peaks for 
which the first upcrossing occurs on the filtering scale Rs is 



o3/2 f f fO fSc-S'ARsJ: 

nvci'^c,Rs)ARs = d^C d^V d5' — |detC| (5z)(r?) OniXs) Pi{w) 

J J J — oo J 5q 

= Jd'^C jdS jdi^ dfi-^^np^{y)Pi{w)ARs , (5.1) 

where w is the 11-dimensional vector of variables w = (/i, v^rj^C,) = (^u, y). Therefore, we can 
write the excursion set peaks multiplicity function as 

hs?{^c,Rs) = ^ nucivc, Rs)'^ = ^ [d^^WfieHifJ-)npi^{y)Pi{w) , (5.2) 

p dvc J^c% J 

where it is understood that -Rs is allowed to vary while the density threshold 5c is kept fixed. 
We can thus think of excursion set peaks as arising from the discrete number density 

nESp{w) = — ^enifi) ^^pk(y) ■ (5.3) 

1/(Tq 

Our calculation is valid for any smoothing kernel. In the special case of Gaussian filtering, 
fi = Rs(J2U and ctq = —RsCrf/ao, and we recover the prefactor oiu/{'yiv) obtained by [24] [the 
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step function 0_f/(/^) = Oh{u) then becomes redundant with ^//(As)]. For practical purposes, 
the excursion set peaks multiphcity function can be computed from 

hsA'^c, Rs) - -y^ ^ , (5.4) 

where Gn"^ is generahzed to 



oo roo 



G^n\imuti,T^c) = df^ifiAfifi) duu''f{u,a)Af{u\uc,ij), (5.5) 
Jo Jo 

with 7tj^ given shortly. Note that the multiplicative factor of f{u,a) remains unchanged 
because the variable is uncorrelated with {v, u, fi). 

Following [23], rotational invariance implies that the 1-point probability density Pi{w) 
be written as 

Pi{w)d^^w = M{v,u,n)diydudn x xi(377^)o?(3r?2) x xl{^Q^)d{K'^) , (5.6) 

where Af{u,u,iJ,) is a trivariate normal distribution and xti^) ^ x^"distribution. The 
cross-correlations between the variables ly, u and ^ are 

{uu) = 71 , {v^) = -ctq = j„f, , (un) = - —a[ = -fu,i , (5.7) 

whereas (z/^) = (u^) = 1 and (/.i^) = ((5')^) = ^o- Owing to the new scalar variable fi, the 
effective local bias relation (3.18) must be generalized to 

Spk{x) = aobiooi^{x) + a2boiou{x) + 6ooi^(a;) 

+ ^cro^200i^^(a;) + (Tocr2biiou{x)u{x) + ^albo2ou^{x) + ajxioV^i^) + o-2XoiC^(a;) 
+ ^boo2lJ''^{x) + aobioiu{x)fi{x) + cr26oiiu(a;)^(a:) H (5.8) 

The bias parameters bijk are defined as ensemble averages over the trivariate Hermite poly- 
nomials constructed from M{u, u, fi), 

<ylcr{hijk = ld^^wnESp{w)Hijk{i',u,fi)Pi{w), (5.9) 

JT-ESP J 

where nESP = /esp/^ is the number density of excursion set peaks. There is no factor of Aq 
in the left-hand side since unlike and u, is not normalized to have unit variance. 

The excursion set peaks multiplicity function scales as /esp exp(— z^^/2)/z/c, which 
should be compared to /(fc) = exp(— f^/2) in the Press-Schechter formalism. The presence 
of a multiplicative factor of z/~^ guarantees that the Gaussian bias factors fefcoo; which involve 
derivatives of Pi with respect to the peak height, agree exactly with those obtained from a 
peak-background split applied to the mass function nh(M) oc I'cfESPi'^c, Rs) of dark matter 
halos. 
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5.2 Non- Gaussian bias 



In light of the results derived in Section 3, a quick calculation shows that the non-Gaussian 
correction to the power spectrum of excursion set peaks is 



'2i '2: 2 2 



+ A^6oo2 - {(yiy hoi - (o-a+i)'6oii 



(5.10) 



The symbol is similar to the spectral moment o"^ defined above, except that the density 
field 5 is replaced by its derivative /u with respect to the filtering scale: 



1 /"OO 



(5.11) 



We shall now proceed analogously to Section 4 and compare the amplitude of the non- 
Gaussian correction (5.10) with logarithmic derivatives of the halo mass function, 



"h(M) = ^-c/ESP(-c,i2.)^y^ 



(5.12) 



constructed from the multiplicity function of excursion set peaks. The subtlety resides in 
handling the derivatives with respect to the smoothing radius Rg. 

5.3 Consistency with peak-background split 

In the presence of a long- wavelength background perturbation, the power spectrum Ps{q) 
of the smoothed density field is rescaled according to (4.14) so that, in the large scale limit 
A: <C 1 considered in this work, the spectral moments cjj and a[ transform as 



o-j I IH 



2 2 

l + e^ + 2e^ 



9 log ag+i 
d log ai 



(5.13) 



where e = 2f^i^{kp)M^~^{k). Hence, both cjj and a'^ equally contribute a term e((TQ+i/(Tj)^, 
whereas a'^ solely induces a correction proportional to 9 log log . Therefore, to com- 

pute the correction proportional to, e.g. ((Tq/cto)^ from peak-background split, we must add 
the logarithmic derivatives of nh(M) with respect to both (Tq and iTq. On reformulating the 
halo mass function in the more convenient form 



-l^cfESp{l^c,Rs) — 



dR 
dM 



(5.14) 
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and taking into account the dependence of the excursion set peaks multiplicity function on 
o"o, we find 



d log nil ^ dloguh 
d log gq d log CTq 



d log /esp d log /esp 2^ + ^ ^^^^ 



d log Pi 
d\og V 



Slog 71 

d log /esp 
log 71 



1 + 



d log 7i.^ 
^ log /esp 
9 log 7^^ 



(5.15) 



where (X) = J Xn-ESpPi/nESP- The last expression involves second derivatives of the 1-point 
probability density Pi{w) with respect to the variables and fi. In particular, 



d log nh ^ S log /esp 
log do Slogo-^ 



+ 1 



^ log /esp 
9 log 7,,^ 



101 



(5.16) 



These relations can now be substituted into the terms proportional to and (cr^)' in the 
square brackets of (5.10), 



'200 



{^D'bioi 



(^0^200 + 27,y^fTo6ioi) ( -f ) + Ti^M^o^ioi 
9 log nil / ^'^o\ 5 log nil / <J<7n 



0-% [d log a a 



2—— 
al V^logcTo 



1 



SlogCTo \ J ^ d\og a'r. \ a'n 



(5.17) 



For fji and o"^, we must consider derivatives of the halo mass function with respect to i?*, 71 
and 7u/i. A straightfoward calculation leads to 



^log^^h ^ d log /esp ^ 9 log /esp 
Slogo-'^ 9 log cr^ 9 log 7^1^ 



'727«/i^011 , 



(5.18) 



and 

9 log nil _|_ 9 log nil 



51og/ESP , 51og/ESP , „91og/ESP o 2/, , N , , 
2^T^ h ^T-; h 2^7-; = 2ai (6110 + Xio j + O-27«M°011 • 



(91og(Ti Slogu'^ (9 log 71 51ogi?- 
Therefore, 

2(7^+1 (6110 + xio) - (fT^+i)'6oii 

= 2(7^ (6110 + xio) + 2o-27«^6oii 
_ d log nil / 5ai \ ^ d log ( ba'^ 



d log 7„ 



(5.19) 



Q + l 



+ O'27tiM^011 



d log fj, 



Q + l 



9 log (Ti 



91og(Ti \ CJi 



+ 



d log fj^ 



0- 



(5.20) 



For (T2, the calculation is somewhat more intricate. On the one hand, taking the derivative 
of the halo mass function with respect to CJ2, we arrive at 



9nii _ 9/esp ^71 ^ 0/esp dRj, ^ 5/esp 97«/^ 



da- 



d-fi da2 dRi, 802 d-fu/i da2 



/esp H [TT 

02 0-2 VKt 



^/2^^ ~^ \ du / da2 



/dPi \ du 



(5.21) 
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The term — 9/esp/o'2 includes the contribution from F(u, w, v) (a factor of —6) and from the 
measure dudvdw (a factor of —3). On the other hand, the terms proportional to (Ja+2 in the 
non-Gaussian bias correction Aci(fc) are 

C^a+2&020 + 2o-^+2X01 = (o-2^020 + 2(72X01 ) • (5-22) 

We can easily convince ourselves that 

2 91og/ESP <91og/ESP o"2 / dPi 



9 log 71 Olog-^uix JESP \ OU 



Therefore, 



2, , o 2 91og/ESP 51og/ESP 0-2 / dPi 

O"2O020 + 20-2X01 ^ 



5 log 71 51og7„^ /esp \ du 

91og/ESP e „o , ^{a), . 

— 8 - 2a„ log ' 71, 7«M, Vc) 



(5.24) 



0=1 



where, in the second line of the right-hand side, we wrote —5 = —8 -|- 3 and used the fact 
that — 9 log /esp /c? log i?* = 3. On multiplying (5.21) by (cr2/nh) and comparing with (5.24), 
we find that 

2 , ,o aiognh /cr2 2\ _ aiognh /5o-2\ 

O"a+2O020 + 2o-a+2X01 = = ^1 • i^'^^) 

91ogCJ2 V / 51ogCr2 V 0-2 y 

Finally, the last contribution to be checked is that from A^. The second-order bias factors 
^002 is proportional to the logarithmic derivative of /esp with respect to Aq, i.e. 

2, _ Slognh 
^°^°°^ = 5biA^- ^^-^'^ 

Therefore, observing that (5Ao/Ao = A^/Aq, the square brackets of (5.10) (i.e. the sum of 
second-order bias coefficients) can be expressed as 

d log fih ^ " 91og?ih / SaA ^ d log nh / ^ d log f 5a2\ _^ d log nh / 5Ao\ 

(5.27) 

where the change in the normalisation amplitude cTg is generally scale-dependent. Note that 
the term ^cto/cto brings down the factor of (Slog aajdXog o-q — 1) found by [14], which leads to 
pronounced effects when the primordial non-Gaussianity is not of the constant-/NL</'^ form 
(see [16]). In the excursion set peaks formalism, this correction appears from the requirement 
that i?s (or M) varies while b{x) = 6c is kept fixed. 



5.4 Non-Gaussian bias and non-universality of the mass function 

We emphasize again that the halo mass function derived from excursion set peaks is not 
universal. In this regards, [15] derived an expression for the non-Gaussian correction to the 
first-order bias which is valid for non-universal multiplicity functions of the form /(ScCTq). 
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In the low-/c limit and for a primordial 3-point function, their result reduces to [omitting a 



factor of 2Af^i^{kp) / Msik)kj 



-2ai 



-la-. 



for clarity] 



d 
dM 



<7. 



-%i^f{6c,cro] 



5 log fJo 



d log (To 



1 + 



dM 



log [^^/(^c, (To)] 



1 dcTQ 

(Tq dM 



. (5.28) 



The main difference with our findings is the presence of a derivative with respect to the halo 
mass rather than the normalisation amplitude erg- Does their formula really differs from 
ours? 

To answer this question, we substitute /esp into the above expression and write the 
derivative with respect to the halo mass as (d/dM) = {d/da(j){dao/dM)+{d/daQ){daQ/dM)+ 
■ ■ ■ . Using (5.13), we eventually arrive at 

^ f log /esp ^ ^ log /esp _ ^ _^ 2 f ^ f^^^ + 1^ ^ f ^log^a _ ^ 
\ aiogcJo Slogo-Q ) al V dlogcj'^ ) erg V^logcro / 

/ ^ log /esp _^ ^ log /esp \ O'a+l ^ 2 f '^^°g/ESp \ / ^ log Ug+l _ \ 

SlogcTi 91ogcr^ / cr? V d\oga[ J erf \ Slogcri / 



' 9 log /esp \ 
, 91ogcr2 J a\ 



+ 



^ 9 log /esp \ A| 
SlogAo ) A2 



(5.29) 



A comparison of the right-hand side with expressions derived above shows that it is exactly 
equal to (5.27) [a rather unsurprising result given that smoothing kernels always appear 
through ((/)]. Therefore, for excursion set peaks, the square bracket of (5.10), (5.27) and 
(5.29) are all equivalent expressions for the amplitude of the non-Gaussian correction to the 
linear halo bias. 

The consistency of the formalism considered here is now established. The non-Gaussian 
correction to the first-order peak bias is a weighted sum over the second-order peak bias 
factors. These nicely combine into a logarithmic derivative of the halo mass function with 
respect to the normalisation amplitude (e.g. dg), as expected from peak-background split. 
The first-crossing condition ensures that the correction found by [14] be present and that 
the ^-independent piece of the Gaussian peak bias factors satisfy the peak-background split 
relation b^oo = (,—^Yn^d!^n\^ld^\. 



6 Discussion of the results 



We shall now draw connections with previous analytical works on the non-Gaussian bias, 
focusing on thresholded regions, and discuss why the standard Lagrangian local bias model 
fails at predicting the correct non-Gaussian bias. 

6.1 Non-Gaussian bias of thresholded regions 

Matsubara [20] worked out non-Gaussian bias corrections for rather generic Lagangian bias 
relations. On large scales, the leading order contribution to the non-Gaussian bias was found 
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to be 



A6i(A;) 



Q2{k) 



1 



C2{q,k- q)Bsiq,k,\k- q\) . 



(6.1) 



2Psik) 2Ps{k) J (27r)3 

In addition, the continuous number density n[x, M) of halos of mass M was defined in such a 
way that the actual halo mass function is recovered under spatial averaging. In the particular 
case of regions with overdensity equal to the critical thresholded 5c (what we refer to hereafter 
as thresholded regions), one has 



n{x, M) 



2p d 
M'dM 



(6.2) 



It can be easily checked that {6H[Ssix) — 6c]) = P{M,6c), the probability that 6s{x) exceeds 
the critical threshold for collapse. To extend the scope of his calculation, Matsubara assumed 
that n{x,M) may be different from dM^ni^s — 6c)- For universal multiplicity functions, the 
second-order renormalized bias parameter was shown to take the form 



C2(fci, fca) = h2{M)W{kiR,)W{k2Rs) + 



l + 5cbi{M) d 



Si 



d log do 



W{kiRs)W{k2Rs) , (6.3) 



where bi and 62 are the usual Gaussian peak-background split bias parameters. Note that 
factors of 9/9 log M appear at all orders. As we will see shortly, they are equivalent to the 
fj,"'{x) present in the effective local bias expansion (5.8). As a result, the function Q2{k) is 
represented by 

1 + (A:) 



Q2{k) = 62X2 (A:) + 



52 9 log (To ' 



where 



Mk) 



<fq 

(2^)3 



W{qRs)W{\k- q\Rs)Bs{q,k,\k-q\) 



(6.4) 



.5) 



Inserting this expression into Q2{k) and rearranging the terms, the non-Gaussian halo bias 
in the large-scale limit becomes 



A6i(A;) « 2fni^M{k)-^a^ 



{5% + 25cbi + 2) + 25-2 {l + 5cbi 



d log a a 
9 log do 



. (6.6) 



In the case of a Press-Schechter multiplicity function, we have the additional simplifications 
5lb2 + 25cbi + 2 = {5l/al)bi and l + (5c6i = 1^2, so that the non-Gaussian bias formula derived 

by [14], 



^bi{k) = 2hi^M{k)-^^ 



Scbl 



d log Q-g 

d log fJo 



1 



(6.7) 



is recovered. However, as was noted in [20], this does not occur in cases other than the 
Press-Schechter mass function (for the Sheth-Tormen mass function [44] for instance). 

Ref. [21] also explored the clustering of thresholded regions, computing 2-point statistics 
from a local bias expansion formulated in terms of Hermite polynomials (as done in [45] for 
Gaussian initial conditions). They showed that such a local bias expansion leads to a non- 
Gaussian bias consistent with peak-background split expectations. This suggests that it 
should be possible to rephrase the calculations of [20, 21] in the formalism considered in this 
paper. 
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n 



6.2 Connection with excursion set peaks 

To emphasize the connection with the local bias expansion (5.8) formulated for excursion set 
peaks, let us go one step further and perform the derivative of the step function with respect 
to to the smoothing scale Rs- The number density of thresholded regions on that filtering 
scale becomes 

n{x,R,) = -'^^eH[5s{x) - 5,] = f^^^SoHx) - . (6.8) 
For a Gaussian density field, the average number density of thresholded regions thus is 

J M do 

(6.9) 

where we have used Bayes' theorem A^(z/c,/i) = M{^\i'c)M{vc) and the conditional average 
{^\i'c) = lu^^l^c = — ctqI^c- Since duc/dRg = —VcicyQ/ao), the average number density of 
thresholded regions in the infinitesimal range [z^o t'c + dvc\ is 

dB 1 

n{v,) = n{Rs)—^ = -2V-' / dfifi^iu^fi) = 2V-^N{v,) . (6.10) 

Hence, we get the Press-Schechter multiplicity function /ps(t'c) = Vnpsii'c) = 2J\f{vc). Note 
that the trajectories 5{Rs) can equally cross the threshold 5c up or down depending on the 
sign of /i. This lead to the so called "could-in-cloud" problem. In our excursion set peaks 
approach, this issue is taken care of upon requiring /_f > [hence the multiplicative factor of 
Ofji^j) in (5.3)], which turns out to be a very good approximation at large smoothing scales 
(where trajectories are nearly fully correlated, see [43]). 

We can now easily convince ourselves that clustering statistics of regions at the threshold 
5c can be computed from the local series expansion 

5y{x) = hio5{x) + hQifi{x) + hi2o5'^{x) + hii5{x)^i{x) + ^6o2/U^(a3) H (6.11) 

with the understanding that, in analogy with discrete peaks, terms involving zero-lag mo- 
ments must be discarded. The bias factors are computed analogously to those of discrete 
density peaks. Since u and /i are normally distributed random variables, the bias parameters 
hij are bivariate Hermite polynomials averaged over all possible thresholded regions, i.e. 

(Tohj = _^ ldudfin{x, Rs)Hij{u, fi)J\f{u, fi) . (6.12) 
n[Ks) J 

Following [23], the biases bij can be evaluated straightforwardly from the series expansion 
[f{eue2)e^^-ob.+e,b,^^^^~^ ^ ^ ^ alU^(^ , (6.13) 

where /(ei, £2) is the exponential factor in the bivariate normal M(y^ fi) with the replacement 
— )• ei and ^ — )• £2, and 
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The difference with bi^ and 6^ defined in (3.14) arises from the fact that the rms variance of 
/i is not normahzed to unity. The first-order bias factors thus are 

-1 r+°° n. f Kir 



510 



1 








1 



















(6.15) 



for the density 5{x), and 



yoi 



(6.16) 



for its derivative ijl{x) with respect to the filtering scale. To derive these results, we took 
advantage of the fact that (/i^|i^c) = (a'I'^c)^ + (A/i^|i/c), where (A/i^| i/c) = Aq — is the 
variance at a fixed value of = Vc- Consequently, the first-order Fourier space bias factor 
ci{k) of thresholded regions is (adopting the notational convention of [20]) 



ci{k)5{k) = bio5sik) + boi^ik) 



1 dW 

biiM)W + 



6c d log (To 



5{k) , (6.17) 



which agrees with the expression found in [20, 21]. A similar calculation shows that, at second 
order, 620 = {^c ~ 3)/o'o5 ^11 = ~l/('^ofo) ^02 = 0, so that we recover the bias factor 
C2(A;i, ^2), (6.3), in the particular case of thresholded regions [for which b2{M) = (z^^ — 3)/(Tq]. 
This shows that clustering statistics of thresholded regions can be computed exactly like those 
of discrete peaks. 

As shown in [23] and in the present work, the effective local bias expansion vastly 
simplifies the calculations for discrete density peaks. Still, computations are somewhat more 
intricate than for thresholded regions. Therefore, the question arises as to whether it would be 
possible to modify the "localized" number density (6.8) so as to obtain, e.g. a Sheth-Tormen 
multiplicity function and the corresponding bias factors, while simultaneously accounting for 
the correct non-Gaussian bias amplitude. A sensible choice is 

n{x,R,) = ^t^W{u,f,)5DHx) - vc] , (6.18) 

where W(i^, /u) is some weight function to be determined through the calculation. Note 
that it cannot depend on u only, otherwise the bias factors computed from (6.12) would be 
the same as those of thresholded regions. The average number density becomes n(fc) = 
2V~^g{ucW{j^c), where 

3(^c) = r^M^c, I^Wii^Wc) = W(^,) - (aI - a',' + + . . . _ 

(6.19) 
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The second equality assumes that W{vc, n) is Taylor expanded around /i = 0. The first-order 
Gaussian bias factor 610 can be derived either from a peak-background split, 



1 dn 1 

^10 



(Ton dvc ctq 



1 g'i^c) 



+ (6.20) 



or from an evaluation of the ensemble average of b^, 



no 



(6.21) 

A detailed calculation yields oi ^02, which suggests that it may be difficult to ensure the 
equality of both expressions except for very specific choices of W. Furthermore, 602 / 
so that, unlike thresholded regions, the local bias expansion will also involves ^^(a;)-terms 
etc. Therefore, one cannot simply set /ps —5- /st in the local bias expansion (6.11) [or, 
equivalently, consider only terms like {d/dM)(yVi x • • • x Wn) as done in [20]]. We will not 
explore this issue any further here. The bottom line is working out a fully consistent biasing 
scheme based on a fitting formulae for the multiplicity function is not trivial. 

6.3 Why does standard local bias fail? 

We are now in a position to address the bi versus 62 issue mentioned in Section 1. As we 
have seen, the peak-background split ansatz is not absolutely necessary to the calculation 
of clustering statistics of discrete density peaks (and thresholded regions). However, it is 
an essential ingredient for the following reasons: i) it can dramatically simplify the calcu- 
lations (through the effective local bias expansion proposed in [23]), and ii) it provides two 
consistency relations that must be satisfied by any realistic (Lagrangian) description of halo 
clustering. Namely, 

• The /c-independent piece b^ {= h^QQ for the excursion set peaks) of the Gaussian bias 
factor CAr(fci, . . . , k^) is equal to a derivative of the halo mass function. 



• In the low-A; limit, the amplitude of the non-Gaussian correction Aci(A;) to the linear 
bias ci{k) is proportional to (any proxy for the normalisation of Ps is acceptable) 

Ac,(fc) oc . (6.23) 

Even though excursion set peaks and thresholded regions satisfy both equalities, it is generally 
difficult to satisfy them simultaneously. For example, changing the number density (6.8) into 
(6.18) in an attempt to improve the agreement with N-body multiplicity functions appears 
to violate both conditions. This suggests that only Lagrangian bias schemes defined through 
a set of constraints imposed on the linear density field can satisfy these peak-background split 
relations. 

Therefore, the fundamental reason why the standard Lagrangian local bias model 5-^{x) = 
bi5{x) + b2S'^{x)/2 + ■ ■ ■ fails at reproducing the correct non-Gaussian amplitude originates 
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from the fact that, unhke the discrete density peaks and the thresholded regions considered 
here, the local bias expansion is not self-consistently computed from a constrained subset 
of the linear density field [this statement holds true also for the approach of [17] which, for 
/nl = 0, reduces to standard local bias]. The constraints must involve 5{x) solely for a local 
bias expansion Sh{x) = J^[6{x)] to give results consistent with the peak-background split 
expectations. This is a very restrictive condition given that, e.g. terms involving powers of 
fj,{x) will appear in the series expansion as soon as any barrier crossing condition is imposed. 
Since the Gaussian bias 61 = (f^ — l)/5c of thresholded regions does not agree with measure- 
ments of the linear halo bias for realistic values of Vc, we are led to the inescapable conclusion 
that a fully consistent modeling of halo clustering statistics with Gaussian and non-Gaussian 
initial conditions requires variables in addition to S{x) and ^{x). In the case of excursion set 
peaks, clustering statistics involve also the curvature of the density field u{x) etc. (note that 
it should possible to include such dependencies in the approach of [19]). The multiplicity 
function /espI^^c) oc exp(— f^/2)/z/c guarantees that the peak-background split consistency 
relations hold for both the Gaussian bias parameters and the non-Gaussian corrections. 

7 Conclusions 

We have computed non-Gaussian corrections arising from the primordial bispectrum to the 
2-point correlation function of discrete density peaks. We have shown that the local bias 
expansion of [23] gives rise to the same result as a direct calculation from the Edgeworth 
expansion. Furthermore, we have generalized the former to the excursion set peaks of [24]. 
In all cases, we have checked that the non-Gaussian correction to the linear peak bias agrees 
with peak-background split expectations. More precisely, the effective local bias expansion 
predicts a non-Gaussian amplitude which is a sum over quadratic bias parameters. All these 
terms nicely combine into a derivative of log with respect to the normalisation amplitude 
CTg (or, equivalently, the halo mass M). This is very encouraging since excursion set peaks 
have recently been shown to reproduce very well the Gaussian mass function and bias of 
dark matter halos [25] . Even though we have focused on the contribution of the primordial 
bispectrum, our analysis can be straightforwardly extended to generic non-Gaussian initial 
conditions, and provide the basis for a computation of non-Gaussian corrections to the halo 
bispectrum. 

Our results shed new light on the widespread local bias expansions and their ability 
to model Gaussian and non-Gaussian halo clustering statistics. We argue that generic La- 
grangian bias expansions cannot properly model these statistics (and, thereby, satisfy the 
peak-background split relations mentioned above) unless the series expansion consistently 
arises from a set of constraints applied to the linear density field (as is the case for peaks or 
thresholded regions). Since the trivial constraint 5s{x) > 6c does not furnish a good fit to 
halo statistics, we are inescapably led to consider additional variables beyond the density 6s 
and, therefore, local bias relations beyond the widespread series expansion 6yi{x) = J^[6s{x)]. 
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